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We investigate propagation of slow-light solitons in atomic media described by the nonlinear A- 
model. Under a physical assumption, appropriate to the slow light propagation, we reduce the 
A-scheme to a simplified nonlinear model, which is also relevant to 2D dilatonic gravity. Exact 
solutions describing various regimes of stopping slow-light solitons can then be readily derived. 
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Recent progress in experimental techniques for the co- 
herent control of light-matter interaction opens many op- 
portunities for interesting practical applications. The ex- 
periments are carried out on various types of materials 
such as cold sodium atoms 0, rubidium atom vapors 
[J solids 0, ll| , and photonic crystals . These 
experiments are based on the control over the absorption 
properties of the medium and study slow light and su- 
perluminal light effects. The control can be realized in 
the regime of elcctromagnetically induced transparency 
(EIT), by the coherent population oscillations or other 
induced transparency techniques. The use of each differ- 
ent material brings specific advantages important for the 
practical realization of the effects. For instance, the cold 
atoms have negligible Doppler broadening and small col- 
lision rates, which increases ground-state coherence time. 
The experiments on rubidium vapors are carried at room 
temperatures and this does not require application of 
complicated cooling methods. The solids are a strong 
candidate for realization of long-living optical memory. 
Photonic crystals provide a broad range of paths to guide 
and manipulate slow light. The interest in the physics 
of light propagation in atomic vapors and Bose-Einstein 
condensates (BEC) is strongly motivated by the success 
of research on storage and retrieval of optical information 
in these media d S S 13, E3, 

Even though the linear approach to describing these ef- 
fects based on the theory of electromagnetically induced 
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transparency (EIT) 12] is developed in detail mod- 
ern experiments require more complete nonlinear descrip- 
tions jl 1| . The linear theory of EIT assumes the probe 
field to be much weaker than the control field. To al- 
low significant changes in the initial atomic state due 
to interaction with the optical pulse, in our considera- 
tion we go beyond the limits of linear theory. In the 
adiabatic regime, when the fields change in time very 
slowly, approximate analytical solutions |l4. Il5l | and self- 
consistent solutions [lfj were found and later applied in 
the study of processes of storage and retrieval [13 . Dif- 
ferent EIT and self-induced transparency solitons in non- 
linear regime were classified and numerically studied for 
their stability 0. As was demonstrated by Dutton and 
coauthors |19| strong nonlinearity can result in interest- 
ing new phenomena. Recent experiments and numerical 
studies 1 1 . 20 have shown that the adiabatic condition 
can be relaxed, allowing for much more efficient control 
over the storage and retrieval of optical information. 

In this paper we study the interaction of light with a 
gaseous active medium whose working energy levels are 
well approximated by the A-schcmc. Our theoretical 
model is a very close prototype for a gas of sodium atoms, 
whose interaction with the light is approximated by the 
structure of levels of the A-type. The structure of lev- 
els is given in Fig. ^ where two hyperfine sub-levels of 
sodium state 3 2 Si/ 2 with F = 1,F = 2 are associated 
with 1 2) and |1) states, respectively Q. The excited 
state j 3) corresponds to the hyperfine sub-level of the 
term 3 2 P 3 / 2 with F — 2. We consider the case when the 
atoms are cooled down to microkelvin temperatures in 
order to suppress the Doppler shift and increase the co- 
herence life-time for the ground levels. The atomic coher- 
ence life-time in sodium atoms at a temperature 0.9/iK is 
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FIG. 1: The A-scheme for working energy levels of sodium 
atoms. The parameters of the scheme are: Wi2/(27r) = 
1772MHz, w/(2tt) = 5.1 ■ 10 14 Hz (A = 589nm), and A is the 
variable detuning from the resonance. 

of the order of 0.9 ms Q. Typically, in the experiments 
the pulses have a length of microseconds, which is much 
shorter than the coherence life-time and longer than the 
optical relaxation time of 16.3ns. 

The gas cell is illuminated by two circularly polarized op- 
tical beams co-propagating in the z-direction. One beam, 
denoted as channel a, is a c^-polarized field, and the 
other, denoted as b, is a cr + -polarized field. The corre- 
sponding fields are presented within the slow-light vary- 
ing amplitude and phase approximation (SVEPA) as 

E = e a g ae i(kaZ-Uat) + g b feifaz-urf) + C ( , 

Here, k a ,b are the wave numbers, while the vectors e , e& 
describe polarizations of the fields. It is convenient to 
introduce two corresponding Rabi frequencies: 

2n a E a 2{i b £ b 

— h — ' = — h — ' 

where [i a ,b are dipole moments of quantum transitions in 
the channels a and b. 

Within the SVEPA, and in the sharp line limit case 
(A = 0), the wave equations for two Rabi-frequencies 
are reduced to the first order PDEs: 

d ( Q a = iv ip 3 ipl, = ii/ ip3ip2- (3) 

Here £ = z/c, t = t — z/c, and vq is a coupling constant, 
which depends on the density of atoms. 
The Schrodinger equation for the amplitudes ipi,2,3 of 
atomic wave function reads 

i 

d T ^ 2 = ~o^ 3 ; (4) 

1 



Here 7 describes the relaxation rate, and we set h = 1. 
We can now exclude the amplitudes of the lower levels 
^1,2 and rewrite Eqs.®, in the form: 

^ T i9 c n fc = |o 6; (5) 

9 T |V3| 2 = -7l^3| 2 -^9f(|fia| 2 + |0 6 | 2 ). 
9 rV3 = - 2^ |V- 3 | ^ (I^| 2 9 <<Pa + W O^ifb). 

Here, (p a ,b,3 are the phases of the fields fl a .b and -03, re- 
spectively. For simplicity, and without a loss of general- 
ity, we assume in Eqs.JSJ that the fields £l a ,b are real, i.e 
fa.b = 0. Therefore, we can choose 933 = 0. Notice that 
the first two equations are wave equations for the fields 
in curvilinear space described by the metric depending 
on the amplitude of the excited state 1P3. 
To make parameters dimensionless, we measure the time 
in units of the optical pulse length t p — Ifis typical for 
the experiments on the slow- light phenomena [l| . There- 
fore, the Rabi frequencies are normalized to MHz. The 
spatial coordinate will be normalized to the spatial length 
of the pulse slowed down in the medium to several me- 
ters per second. According to the linear theory, the 

n 2 

group velocity of slow-light pulse is v g « Here 
f^o is a magnitude of the controlling field required in 
EIT experiments. Typically, this field has a magnitude 
of order of few megahertz. So, we choose f2o = 3 and 
v g = 10 _7 c as representative values reported in experi- 
ments with BEC of sodium atoms. Hence the pulse spa- 
tial length is l p = v g t p = 30/im, and £ is normalized to 

10~ 13 s. In the dimensionless units, the coupling constant 

n 2 

^0 = -f = 4.5. 

In the absence of relaxation, i.e. for 7 = 0, the system of 
equations Eqs.©,® is exactly solvable in the framework 
of the inverse scattering method (IS) [2II2I I2I |24|. 
In the present work we provide an elementary method 
to derive slow-light solitons for the case of an arbitrary 
controlling field. 

In the context of slow light phenomena, the system is 
assumed to be initially in the following stationary state: 

fi o =0, n b = fi(r), hu = |i>. (6) 

Notice that the state |1) is a dark-state for the control- 
ling field f2(r). This means that the atoms do not in- 
teract with the field O(t) created by the auxiliary laser. 
The configuration Eq.© above corresponds to a typical 
experimental setup (see e.g. 0,00)- 
We intend to study the dynamics of coupled atom-field 
modulations in the A-type model, which can preserve 
their spatial shape to a large extent while propagating 
in the media. We consider such solutions as a gener- 
alization of the dark-state polariton [2^|. In the linear 
theory the probe field only appears in f2 a , whereas in the 
nonlinear theory it also forms an inseparable nonlinear 
superposition with the controlling field in the channel b. 
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However, in both cases the Rabi-frequency f2 a describes 
probe field modulations. Indeed, the field in the channel 
a induces atomic transitions from the state |1) to the ex- 
cited state 1 3). On the other hand, the amplitude of the 
excited state drives the field £l a ■ From the results of lin- 
ear and nonlinear theories of electromagnetically induced 
transparency (EIT) [U d Q El HI HI we can infer 
a physically plausible assumption that the population of 
the upper level is proportional to the intensity of the field 
in the probe channel a, i.e. |£! a | 2 ~ iV'sl 2 - I n the present 
work we assume that this observation is relevant for slow 
light phenomena. Therefore, we postulate that 



l^3| 2 



-m a \\ 



(7) 



where k is an arbitrary parameter. 

We emphasize that the imposed constraint Eq. Q reduces 
Eqs. JSJ) to a simplified nonlinear system, which provides 
adequate descriptions of the slow light propagation. In 
this sense the relation Eq.0 is central for the present 
work. As we show below this condition is sufficient and 
necessary for the slow-light solitons to exist. Introducing 
new notations: |fi a | = e~ p ,£!{, = ?y, we find from Eq.Q 
together with Eq. Q that the field p satisfies the Liouville 
equation 



d <T p = -ke- 2 ?. 



together with the constraint 



d( T 7] + d T p d^rj = k e 2p i], 



and an auxiliary equation 



(4k(d T + 7 ) + »c) e ~ 2p + d cn 2 = °- 



(8) 



(9) 



(10) 



It is interesting that the Liouville equation Eq.JHJ) ap- 
pears in 2D gravity j2?J and describes a gravitational 
field defined by the metric g a b = e~ 2p ^ a b, where 7 a b is 
the 2-dimensional Minkowski metric. This connection to 
2D gravity is further emphasized by the observation that 
any solution of dilatonic equations [2?! satisfies the sys- 
tem of slow-light Eqs.JSJ), J^J. The dilatonic equations 
have the following form 

d v V{r]) _ 2 V{rj) _ 2 

d CT p '— — e p = 0, d Cr r] + —— e p = 0, (11) 



together with the constraint 

d (T r) + 2 d T p d ( r) = 2kA e~ 2p 



(12) 



Here A((, r) is an arbitrary source term and V(rj) — 
4fc (^l — 77) plays the role of the dilatonic potential. For 
the realization ^4(C, t) = 9 r m(r) with an arbitrary func- 
tion m(r), the equations Eas. ((TT|l . (fT2^l . and Eq. lfTU|) for 
7 = can be readily solved, viz. 



a c A + (pa T a (t) 

(l-fcA + A_) 2 



p = -i log 

A +(0 = -iexp[-8e feC], 
A_(t) =exp heo J 



(It 



2 (r) + l 

t] = 2(d T m — md T p). 



The original fields Q a 5 then read: 



n a = 



2e 



=sech(^), 



= ^Stanh(^) + i^^. 



(13) 
(14) 
(15) 
(16) 

(17) 

(18) 
(19) 



where e is a real arbitrary constant defining the am- 
plitude of slow- light soliton. The background field f2(r) 
reads: 



O(r) 



i9 T m(r) — 2e m(r) 
m 2 (r) + 1 



(20) 



The function f2(r) describes the controlling field, which 
governs the dynamics of the system. The time depen- 
dence of this function is determined by modulation of 
the intensity of the auxiliary laser. As can be readily 
seen, the velocity of the slow-light soliton reads 



1 



1 



4fc m 2 (t) + 1 



(21) 



For a constant controlling field fi(r) = £lo, and in the 

n 2 

simplifying approximation <sC 1, the group velocity 

of the slow-light soliton conform to the result of linear 
theory: 



, "0 
'2vq' 



(22) 



Expression Eq. i12_?l) immediately suggests that the sig- 
nal stops, when Qq = 0. Therefore, this expression is 
the main motivational source for the works on slow-light 
solitons (see H^,|2^ and references therein). We envis- 
age the following dynamics scenario. We assume that 
the slow-light soliton was created in the system before 
the time r = and is propagating on the background 
of the constant controlling field f^o- Suppose that at the 
moment r = the laser source of the controlling field is 
switched off. We assume that after this moment the back- 
ground field will decay reasonably rapidly, as described 
by a "switch-off" function /(r). The front of the van- 
ishing controlling field, described by the function /(t), 
will then propagate into the medium, starting from the 
point C — 0, where the laser is placed. The state of the 
quantum system Eq.© is dark for the controlling field. 
Therefore the medium is transparent for the spreading 
front of the vanishing field, which then propagates with 
the speed of light, eventually overtaking the slow- light 
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soliton and stopping it. To realize this scenario, we as- 
sume the controlling field f2(r) to be constant O for 
negative r and a r dependent switch-off function f(r) 
for positive t, i.e. f2(r) = ilo<d(— r) + /(t)O(t). Here 
O(t) is the Heaviside step function, while /(0) = Slo • In 
this setting, the distance that the soliton travels until full 
stop is 



C = 



dr 



Ak J m 2 (r) + 1' 



(23) 



This distance designates a geometrical point in the 
medium, where the slow-light soliton disappears, writ- 
ing itself into the medium as a standing memory bit in 
the form of a localized polarization cluster. 
From Ea. (|20|l a number of exactly solvable regimes for 
the stopping of the slow-light soliton can be identified. 
For m(r) = e aT and r > 0, a > 0, we obtain /(r) = 
(f — £ o) sech(ar) and 

C= -!-ln2 . 

8ak 



Hence, the slower the field decays to zero, i.e. for smaller 
a, the longer the distance that the soliton travels in the 
medium is. Another physically interesting case /(r) = 



n 



was discussed in |2£ 



Discussion. In this paper we derived the slow-light soli- 
ton from the sufficient condition Eq.(0. In fact, the in- 
verse scattering analysis as applied in |2y| to Eqs.©,® 
shows that for slow-light solitons a stronger condition 
holds, namely 



VA3 = - 



2 A 



(24) 



where A is a complex number parameterizing the soliton 
(in our case A = ieo). This means that the condition 
Eq.Q is also the necessary condition for the slow- light 
solitons to exist with k = 8 | A ^' A p . 

In a forthcoming publication we will explain a fascinating 
analogy between the stopping of a slow-light soliton and 
the formation of a black hole in 2D gravity. 
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